The pressure dependence of electron-phonon coupling in the organic superconductor 
K-(BEDT-TTF)2Cu(SCN)2: A comparison of high pressure infrared reflectivity and 

Raman scattering experiments. 
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We determine the pressure dependence of the electron-phonon coupUng constant in k-(BEDT- 
TTF)2Cu(SCN)2 by comparison of high pressure Raman scattering and high pressure infrared re- 
flectivity (IR) measurements. We use comparison of the IR reflection spectrum from deuterated 
and protonated samples to aid the deconvolution of several overlapping phonon modes. These cou- 
pled modes are modelled with a Green's function to extract the linear pressure dependence of the 
individual modes. The Raman active molecular vibrations of the BEDT-TTF dimers stiffen by 
0.1-1 %GPa~^. In contrast, the corresponding modes in the IR spectrum are observed at lower 
frequency, with a pressure dependence of 0.5-5.5 %GPa~^, due to the influence of the electron- 
phonon interaction. The stronger pressure dependence of the central C=C mode of the BEDT-TTF 
molecule in the IR, is discussed. Our analysis suggests that reduction of electron-phonon coupling 
under pressure does not account for observed suppression of superconductivity under pressure. 



K-(BEDT-TTF)2Cu(SCN)2 is one of the best char- 
acterized organic superconductors [0. It is a highly 
anisotropic material with a quasi-two dimensional band 
structure, whose Fermi-surface topology has been deter- 
mined by magnetotransport experiments 0, §|. At am- 
bient pressure k-(BEDT-TTF)2Cu(SCN)2 is a supercon- 
ductor with a transition temperature of ~ 10.4 K. 
Tc decreases upon the application of pressure until, at 
pressures P exceeding 0.5 GPa, superconductivity is fully 
suppressed ||^, ^. The effective mass, m*, derived from 
magnetic quantum oscillation measurements, decreases 
linearly with pressure up to 0.5 GPa; above this pressure 
the magnitude of dm*/dP is strongly reduced [Q. In 
contrast, the effective mass derived from optical measure- 
ments, ruopt, decreases approximately linearly through- 
out this pressure range Q . The coincidence of a "kink" in 
the pressure dependence of m* with the pressure above 
which superconductivity is suppressed and the absence 
of a "kink" in the pressure dependence of mopt, indicates 
that the interactions parameterised by m* are connected 
to the superconductivity. 

The key question to address is the effect of these in- 
teractions, i.e. what is the dominant pairing mechanism 
for superconductivity in this material? In this paper we 
compare infrared (IR) Q and Raman scattering ||] mea- 
surements under pressure to determine the role of the 
electron-phonon interaction. This is possible because the 
IR measurement probes the molecular vibrations dressed 
by the electron-phonon interaction ^, |^, whereas non- 
resonant Raman measurements probe the bare mode fre- 
quencies [0, H]. Our recent high-pressure IR study 
of this material revealed anomalously large mode stiffen- 
ing for the molecular vibrations most strongly coupled 
to the electronic excitations. However, determination of 
the pressure dependence of the most strongly coupled 
mode, at around 1300 cm~^ in the IR spectrum, is not 
straightforward due to the mixing of several modes with 
varying pressure dependence. In the current paper we 
use the comparison of ambient pressure reflectance data 



from protonated and deuterated samples to establish a 
coupled oscillator model. This model is then employed to 
determine pressure induced frequency shifts of the inter- 
acting modes. By employing a dimer charge-oscillation 
model 0, |[ |§1 , comparison of the pressure induced fre- 
quency shifts observed in both the IR |4| and the Ra- 
man B spectra determines the pressure dependence of 
the electron-phonon interaction. Four modes are observ- 
able in both the high-pressure Raman and IR spectra 
(see Table ^ . They are labeled with subscripts indicat- 
ing the atoms/bond predominantly involved in the vi- 
bration pO) . In order of increasing frequency they are: 
the C-S mode originating from the BEDT-TTF 60 B3g 
asymmetric vibration, the central C=C mode originating 
from the BEDT-TTF 3 Ag symmetric vibration and two 
Cu(SCN)2 anion modes. 

Owing to its strong electron-phonon coupling |Q, we 
are particularly interested in extracting the pressure in- 
duced frequency shift of the central C=C mode. Dimer- 
ization in the K-phase makes this mode IR active ^ . 
Electron-phonon coupling softens this mode from its 
bare, Raman active, frequency of 1470 cm~^ to around 
1290 cm~^ along the 6-axis and around 1220 cm~^ along 
the c-axis. This brings the mode into resonance with the 
end carbon-hydrogen (C-H) oscillations of the BEDT- 
TTF molecule @. Upon deuteration the C-H(D) modes 
are softened sufficiently to reveal the true line shape of 
the central C=C mode. Comparison of the IR spectrum 
from the protonated and deuterated salts illustrates the 
structure of the C=C mode due to this coupling (see 
Fig. |l|b and ||d); dips occur where the stronger C=C 
mode loses spectral weight to the weaker C-H modes it 
is driving. 

Although mode mixing should be treated from a quan- 
tum mechanical point of view, it has been shown [p^ that 
a classical approach, such as described below, can yield 
information regarding the individual mode frequencies. 
From the equations of motion (Fourier components) for 
a set of N coupled oscillators [Q, a Green's-function 
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matrix approach [|15| is employed to model the dielectric 
response function 

e{io) = Q,{G{Lu)}-'n^, (1) 

where flp is a vector containing each oscillator's spectral 
weight, f2j is its transpose and G is a square matrix of 
dimension TV, with TV the number of oscillators. The 
diagonal terms in G are of the form 

Gaa{uj) = iujl+tiUra-LU^), (2) 

specifying each Lorentzian oscillator of frequency uja and 
damping r^. The off-diagonal terms in G provide the 
coupling and are of the form 

Gabiuj) = {Al, + iu;Tab). (3) 

This coupled oscillator model is employed to analyze our 
room temperature high-pressure IR reflectance [Q. The 
aim is to extract the pressure-induced frequency shift 
of individual modes from the apparently non-linear shift 
(predominantly in the c-axis response) of the reflectance 
peaks as a function of pressure (the crosses in Fig. [l^ and 
Fig. 0c). 

To restrict the possible parameter space encountered, 
a simple (uncoupled) Drude-Lorentz oscillator fit to the 
ambient-pressure room-temperature IR spectrum of the 
deuterated salt was used as a starting point, (the dot- 
ted line in figure ^ and |l|d) . In the absence of coupling 
the asymmetric line shape of the central C=C mode for 
both b and c axes is most accurately reproduced when 
the Drude term is heavily suppressed, i.e. damped to 
the point that it hardly contributes to the reflectance 
in this frequency range. The terms describing the back- 
ground and C=C mode were then held fixed and free 
parameters describing the two C-H modes and their cou- 
pling terms were introduced. A least squares fitting pro- 
cedure was employed, yielding the parameters necessary 
to model the interacting modes in the protonated salt at 
room temperature and ambient pressure (the dashed and 
solid line in figures |^b and |]d respectively) . 

For both b- and c-axis polarizations and both the 
protonated and deuterated samples the damping and 
strength of the two C-H(D) modes tended to the same 
values, ^H(D) ^ 16 cm^^ and i^p^H{D) ^ 0, two 
sharp modes with no infrared activity of their own. The 
couphng strengths, A (see equ.(||)), are hsted in table |, 
with the cross damping term Tat also tending to zero 
for all modes. With all parameters except the three 



Parameter b- 


axis value 


c-axis value 








AcccHi 


260 


230 


AcC,CH2 


260 


230 


AcHi,CH2 


160 






TABLE I: mode coupling parameters. 
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FIG. 1: (a) and (c) diamond/sample reflection surfaces generated 
from the coupled oscillator model for the h and c axis respectively; the 
reflectance scales are adjacent to the vertical axes of (b) and (d). Super- 
imposed are the frequency of reflectance peaks (crosses) uncoupled 
mode frequencies (dotted lines) and the peaks of the coupled oscillator 
model (dashed lines), (b) and (d) The ambient pressure reflection spec- 
trum for polarization parallel to the b and c axes respectively: solid and 
dashed lines for the protonated salt (measurement and flt), dash-dot 
and dotted lines for the deuterated salt (measurement and flt). 



mode frequencies held constant, the pressure-induced fre- 
quency shifts of the reflectance peaks were fltted. The 
data could be successfully modelled with only first order 
(linear) pressure shifts included (see Fig. |l|a and |l|c and 
table 0). 

Including a pressure-dependent Drude response to ac- 
count for the increase in background reflectance in this 
spectral region has a small effect (« -1-0.5 %) on the ob- 
served peak frequencies in the reflectance model. Its 
main effect is to broaden and distort the line shape of 
the phonon modes towards higher frequency, as experi- 
mentally observed at high pressure Q. 

Within the framework of a dimer charge oscillation 
model 1^ it is possible to calculate a dimensionless 
electron-phonon coupling constant. A, using a mode's IR 
frequency, wir, and the corresponding Raman frequency, 
wr, given knowledge (see Table ||) of the energy, ujqt, of 
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Mode Raman Infrared b-axis 

cm-^ + (%GPa-M cm^i + (%GPa~i' 



Infrared c-axis 
cm-' + (%GPa-') 



CJCT 
WCS 
WCC 

WCDi 
WCD2 
U-'anionl 
<^anion2 



886.2 + 0.85 
1467.7 + 0.4 



2064.6 + 0.1 
2106.3 + 0.2 



2910 + 4.0 
883.5 + 0.71 
1290 + 2.5 
1181 + 0.5 
1290 + 0.5 
1027 t 
1116 t 
2067.4 + 0.1 



2390 + 4.0 
873.6 + 1.0 
1210 + 5.5 
1177 + 0.5 
1281 + 0.5 
1027 t 
1122 1" 
2065.6 + 0.15 
2109.3 + 0.2 



t 



No infrared pressure data arc available for the deutcrated salt. 



TABLE II; Raman and IR frequencies and pressure shifts 
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The larger observed Awg for the anion mode in the 6-axis 
response (Table 0), for which the lattice mode is stiffer, 
is a further confirmation of the lattice mode assignment 
given in ||]. 



Mode 


Xb 


d 111 Af, 

dP , 
%GPa-i 


A. 


dlnAe 

%GPa-i 




0.01(1) 


47.9 


0.02(1) 


-11.1 




0.17(1) 


-11.9 


0.20(1) 


-17.3 


^anioiil 














t^anion2 













TABLE III; dimensionless electron-phonon coupling con- 
stants and their pressure derivatives. 



the coupling charge transfer band ^ 

A- 



-'IR 



CT 



(4) 



'^CT 



(the C=C mode, being an asymmetric combination of Ag 
molecular modes, couples to the intra-dimer charge trans- 
fer [ p^ ) . The dominant source of error in this calculation 
is due to the width of luq^. The pressure derivative of 
equation (El), 



dlnA 
dP 



+2- 



^IR 



^IR 



d In LUji d In cjir 



dP 



dP 



CT 



d In cjcx d In uj^ 



dP 



dP 



(5) 



has a weaker dependence on the value of ujct- However 
we found that the dominant error in (||) arises from the 
difficulty in determining an accurate pressure dependence 
of the broad CT band. Constraining the pressure depen- 
dence of ujcT to 4±4 %GPa~^ for both b- and c-axes 
Q, results in an error in the pressure derivative of the 
coupling constant for the C=C mode of ±3.5 %GPa~^ 
for the 6-axis and ±5 %GPa^^ for the c-axis. The er- 
ror in the pressure derivative of the coupling constant 
for the C-S mode is larger (8 %GPa^^ for the 6-axis and 
28 %GPa~^ for the c-axis) because Acs is so small. The 
Cu(SCN)2 anion modes do not exhibit any evidence of 
electron-phonon coupling, i.e. the difference between the 
infrared and Raman frequencies for the anion modes is 
accounted for by purely vibrational coupling. The zone 
center frequency separation, Awg, between the symmet- 
ric (Raman active) and asymmetric (IR active) combina- 
tions of the molecular vibrations, cj,noi, is determined by 
the frequency of the optical branch of the lattice modes, 

^lattice [H H , 



^lattice 
2 Wmol 



(6) 



Kozlov et al. calculate that the antiphase combi- 
nation of Bag molecular modes couple to charge trans- 
fer perpendicular to the intra-dimer direction. This sug- 
gests that the C-S mode couples to intra- not inter-band 
electronic transitions. It should be noted however, that 
the electron-phonon coupling constants for the C-S mode 
have been calculated using the same wct as the C=C 
mode because it is impossible to distinguish the contri- 
butions to the IR spectrum from inter- and intra-band 
transitions. 

It has been shown previously that the usual 
electron-acoustic-phonon interaction mechanism is un- 
able to account for the magnitude of the electron-phonon 
coupling constant or the large pressure dependence of Tc 
in the BEDT-TTF superconductors. A further refine- 
ment [Q is to include the attractive interaction medi- 
ated by the Ag molecular modes, with the total elec- 
tron phonon coupling constant, Atot, given by a Yamaji 
sum over the individual Ag molecular modes. The en- 
ergy scale for the interaction is still set by the Debye 
frequency, 9 [|8[ |l^ . Caulfield et al. |^ have previously 
determined 8 « 40cm^^ by fitting the effective mass 
dependence of Tc with a linearised Eliashberg equation 
using an Einstein density of phonon states, 5{Q). High 
and low temperature specific heat measurements [E0| |2l[] 
yield values of ranging from 38cm~^ to 140cm~ . 

Calculations of Atot @ give values ranging from 0.3- 
0.45 [^2[ Based on the vibrations we sampled, we 
assume Atot to have a pressure dependence similar to 
that of the strongly coupled C=C mode, i.e. of the or- 
der of -17 %GPa-i, with an upper limit of -20%GPa-i 
used in this calculation. The weak-coupling BCS formula 
gives a satisfactory explanation of the ambient pres- 
sure Tc and accurately describes the effective mass 
dependence of the superconducting transition tempera- 
ture |, ||. 

The pressure derivative of the weak-coupling BCS for- 
mula provides a convenient parameterization of '"p^*" in 
terms of 



din Tc dine 



dP 



dP 



(A~/x*)2 



dA d^i* 
dP ^ dP 



(7) 
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With = [ln(^)]2, e « 90±50 cm-i |, || 

and using the average pressure induced stiffening of the 
Raman active lattice modes « +13 %GPa~^ |5| for 
the only unknown is the pressure dependence of 

the Coulomb pscudopotential, ^p-. Calculations of /^* 
1^, indicate a small positive value, a possible indica- 
tion that direct interactions between the quasiparticles 
are involved in the pairing. There is also sufficient uncer- 
tainty in the parameters necessary to calculate the pres- 
sure dependence of the dimer-site Coulomb repulsion 
that no reliable prediction can be made. We have thus de- 
cided to ignore in the above calculation. This gives 
« -40±32 %GPa-\ which is far from the ex- 
perimentally observed value of ^'^p^" ~ -200 %GPa~^ 
^. As can be seen from (0), 6 directly scales with 
the pressure dependence of Tc- To obtain the observed 
rapid fall of Tc with pressure from only the decrease in 
the electron-phonon coupling constant, requires to be 
of the order of the C=C mode frequency, « 1500 cm~^. 
Such a value is utterly inconsistent with the 10 K super- 
conducting temperature and the unconventional isotope 
shift observed upon carbon substitution p7| . 

In conclusion we have shown that it is possible to 
recreate the line shape of the resonance between the 
central C=C mode and the H (or D) modes in proto- 
nated (or deuterated) k-(BEDT-TTF)2Cu(SCN)2. This 



was achieved by coupling the strongly infrared active 
C=C mode to two C-H modes which possessed negligi- 
ble IR strength of their own. We have also shown that 
the apparent non-linearity of the pressure dependence of 
the modes is due to anti-crossing of the mixed modes. 
We compare high-pressure Raman scattering and IR re- 
flectivity data enabling the pressure dependence of the 
electron-phonon coupling strength to be evaluated for 
modes observed in both spectra. Using the weak cou- 
pling limit of BCS theory we have been able to compare 
the pressure dependence of the electron-phonon coupling 
constant and the pressure dependence of the supercon- 
ducting transition temperature. This casts considerable 
doubt on whether this material is a simple BCS super- 
conductor because the characteristic energy of the pairing 
interaction would have to be of the order of the highest 
frequency molecular modes, a value inconsistent with the 
10 K superconducting transition temperature. This is an 
indication that electron-electron interaction may be play- 
ing a significant role in this material's superconducting 
mechanism. 
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